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ABSTRACT 

The planetary nebula NGC6881 displays in the optical a quadrupolar morphology consisting of 
two pairs of highly collimated bipolar lobes aligned along different directions. An additional bipolar 
ejection is revealed by the hydrogen molecular emission, but its wide hourglass morphology is very 
different from that of the ionized material. To investigate in detail the spatial distribution of molecular 
hydrogen and ionized material within NGC 6881, and to determine the prevalent excitation mechanism 
of the H2 emission, we have obtained new near-IR Br7 and H2 and optical Ha and [N ll] images, as 
well as intermediate resolution JHK spectra. These observations confirm the association of the 
H2 bipolar lobes to NGC 6881 and find that the prevalent excitation mechanism is collisional. The 
detailed morphology and very different collimation degree of the H2 and ionized bipolar lobes of 
NGC 6881 not only imply that multiple bipolar ejections have occurred in this nebula, but also that 
the dominant shaping agent is different for each bipolar ejection: a bipolar stellar wind most likely 
produced the H2 lobes, while highly collimated outflows are carving out the ionized lobes into the 
thick circumstellar envelope. The asymmetry between the southeast and northwest H2 bipolar lobes 
suggests the interaction of the nebula with an inhomogeneous interstellar medium. We find evidence 
that places NGC 6881 in the H ll region Sh 2-109 along the Orion local spiral arm. 
Subject headings: infrared: ISM: continuum — ISM: molecules — planetary nebula: individual (NGC 
6881) 



1. INTRODUCTION 

Asymmetry is common among planetary nebulae 
(PNe) and it comes in a large variety of shapes and 
morphological features. On a large scale, the departure 
from symmetry goes from the mild asymmetry of ellip- 
tical PNe to the strong asymmetry displayed by bipolar 
PNe. Asymmetry in PNe is also revealed by an assort- 
ment of small-scale morphological features including the 
symmetric fast low-ionization emission regions (FLIERs) 
of elliptical PNe, as well as point-symmetric collimated 
outflows and their associated blowouts and bow-shock 
structures. The asymmetry in PNe has been linked 
on many occasions to binarity (e.g., Sokerl [2QQ6I ) . but 
there is no conclusive evidence supporting such connec- 
tion (Sorensen & Pohacco 2004; Bond 2000; Soked ll997l . 
and references therein). On the contrary, the onset of 
bipolarity in PNe seem s to require the evolutio n of a 
massive progenitor (e.g. JCorradi fc Schwarzlll995l ). 

PNe that show asymmetric structures along different 
directions are of great importance, as the change in the 
direction of the ejection can be associated with the pre- 
cession of the progenitor star in a binary system. This 
class of PNe is not homogeneous at all, and includes 
PNe with multiple collimated outflows along different di- 
rections (Corradi et al. 1 997|) or with p recessing jet-like 
features (Mirand a. Guerrero fc TorrellesI 11999). starfish- 
shaped and multipo lar PNe that show point-synimetri c 
bow-shock features (jSahai fc Traugerlll998l : ISahaill2000[ ). 

^ Visiting Astronomer, Instituto de Astrofisica de Andalucia. 



and quadrupolar and polypolar PNe with different sets 
of bipolar lobes aliened along differ ent symmetry axes 
( Manc hado. Stanghellini fc Guerrero! 1 19961 : iLopez et al.] 
1998). Polypolar and quadrupolar PNe are especially 
interesting as multiple hourglass structures imply recur- 
rent bipolar ejections which are particularly difficult to 
interpret in th e framework of the int eracting stella r win d 
model (jKwok. Purton fc Fitzgeraldiri978l : lBalicklll987h . 

One of the most intriguing PNe with multiple bipolar 
lobes is NGC 6881. Originally classified as a quadrupo- 
lar PN, based on optical narrowband images and long-slit 
echelle spectra (Guerrero fc Manchado 1998), NGC 6881 
presents two pairs of highly collimated bipolar lobes with 
very similar, although not coincident, symmetry axes. 
The southwest lobe displays a loop-like feature that is 
highly reminiscent of a precessing collima ted outflow, 
but its measured expansion velocity is low. iKwok fc^ul 
(2005) showed that the dense equatorial ring has recently 
changed its orientation, being aligned with the youngest 
pair of bipolar lobes. 

Therefore, there are many signs suggesting that pre- 
cession and sequential events of bipolar ejections have 
occurred in NGC 6881. The spatial distribution of the 
molecular hydrogen in this nebula, as revealed by near-IR 
narrowband H2 (1-0) S(l) images (|Guerrero et al.|[2QQQl ). 
adds a new twist. Molecular hydrogen emission is de- 
tected in wide hourglass bipolar lobes that extend much 
farther than the ionized bipolar lobes. The II2 bipolar 
lobes may represent a bipolar ejection unrelated to the 
formation of the two pairs of ionized bipolar lobes. 
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TABLE 1 
Narrow-Band Imaging of NGC 6881 



TABLE 2 
NIR Spectroscopic Observations of NGC 6881 
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9 
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To study carefully the spatial distributions of ionized 
material and molecular hydrogen in NGC 6881 and to 
investigate the excitation mechanism of the H2 molecule 
throughout the nebula, we have obtained new H2, Br7, 
[N n], and Ha images, and JHK intermediate-resolution 
long-slit spectroscopic observations. A description of the 
observations is presented in §2 and the results are given 
in §3. The results are discussed in §4, and the conclusions 
and a short summary are presented in §5. 

2. OBSERVATIONS 
2.1. Imaging 

Narrowband near-IR images of NGC 6881 were ob- 
tained during 2006 September 9 using LIRIS (Long-Slit 
Intermediate Resolution Infrared Spectrograph) at the 
Cassegrain focus of the 4.2 m William Herschel Tele- 
scope (WHT) on Roque de Los Muchachos Observatory 
(ORM, La Palma, Spain). The detector was a Ikxlk 
HAWAII array sensitive in the spectral range from 0.8 to 
2.5 /im. The plate scale is 0'.'25 pixel" ^ and the field 
of view (FOV) is 4.'27x4.'27. The narrowband filters 
isolated the H2 (1-0) S(l) 2.1218 //m and Br7 2.1658 
/xm emission lines. An additional narrowband filter (Kc) 
centered at 2.270 jivn was used in order to subtract the 
continuum contribution from the line emission detected 
through the H2 and Br7 filters. The central wavelength 
and bandwidth of these filters are listed in Table 1. 

We obtained series of 100 s exposures on each filter 
for the total integration times given in Tab. 1. The tele- 
scope pointing was shifted by a few pixels between each 
exposure, rastering the nebula to different locations on 
the detector. Each series of observations on the object 
was followed by a similar series of observations on ad- 
jacent blank sky positions. Individual exposures were 
fiat-fielded and dark-corrected, and the sky frames were 
combined to obtain a master sky image. The images 
on the nebula were subsequently sky-subtracted, shifted, 
and averaged. The continuum contribution to the H2 and 
Br7 emission was removed using the continuum Kc im- 
age. Since the H2, Br7, and Kc filters have very similar 
bandwidth, we do not find it necessary to scale the im- 
ages, although differences in the point-spread function 
(PSF) of stars in the field of view required us to de- 
grade the PSF of the H2 image to match these of the 
Br7 and Kc images. The continuum subtracted H2 and 
Br7 images and the Kc image are presented in Figure [H 
The spatial resolution, as determined from the FWHM 
of stars in the FOV, is -0'.'8. 

Narrowband images in the [N ll] and Ha emission lines 
were obtained on 2006 June 30 using ALFOSC (An- 
dalucia Faint Object Spectrograph and Camera) at the 
2.56 m Nordic Optical Telescope (NOT) in the ORM. 
The camera is a 2048x2048 CCD with a plate scale of 



0:'19 pixel-^ and a FOV of 6f5x6.'5. The [N 11] and 
Hof images of NGC 6881 are also presented in Fig. [TJ 
The spatial resolution, as determined from the FWHM 
of stars in the FOV, was 0'.'9. 

2.2. Spectroscopy 

Intermediate-resolution JHK long-slit spectroscopic 
observations were performed on 2004 July 2-4 using the 
NICS (Near Infrared Camera Spectrometer) at the 3.5 
m Telescopio Nazionale Galileo (TNG) at the ORM. 
The NICS is a multimode instrument for IR observations 
(0.9-2.5 /im) that uses a HAWAII 1024 x 1024 array as 
detector. We used the LF camera, providing a plate scale 
of 0'.'25 pixer\ and the 0.75 slit, with a width of 0'.75 
and a length of 4^ Observations were obtained using the 
JH and K^ grisms and the long-slit was placed along the 
central star at P.A.s 113° and 137° (Table 2). In order to 
provide the means for subtracting the sky emission, the 
nebula was placed at different positions along the slit. 

The spectra were reduced using IRAF^ routines of the 
noao.twodspec and noao.onedspec packages. The data 
were fiat-fielded using dome fiats, and the sky contri- 
bution was removed using sky spectra obtained at the 
same time as the nebular spectra. We used telluric sky 
lines for the wavelength calibration. The standard IR 
stars of spec tral type A- AO SAO 15832, SAO 48300, and 
SAO 72320 (|Hunt et al.l ll998) were used for the fiux cal- 
ibration. We obtained the sensitivity function by com- 
paring the spectra of these stars with a blackbody model 
of temperature 9480 K, similar to the effective tempera- 
ture of AO stars. Finally, the telluric absorptions in the 
one-dimensional spectra of NGC 6881 were removed us- 
ing the IRAF task "telluric", applying it to stars with 
featureless spectra in the JHK range. 

3. RESULTS 

3.1. Morphology 

The notable differences in the spatial distributions 
of ionized material an d molecular hydroge n within 
NGC 6881 described bv iGuerrero et all (|200 0l are spe- 
cially highlighted in the H2, Ha, and [N ll] composite 
picture shown in Figure [21 In NGC 6881, we can distin- 
guish the following components. 

3.1.1. The Equatorial Torus at the Central Region 

The central regions of NGC 6881 contain both molec- 
ular and ionized material. This region can be de- 
scrib ed as a clumpy torus-like structure (|Kwok fc Sul 
I2OO5I ) that is expanding at a moderate velocity 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation 
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Fig. 1. — Images of NGC 6881 in the Brj (top, left), Kc (middle, left), and [N ll] (bottom, left) narrowband filters. The images in the 
H2 (top, right), [N ll] (middle, right) and Ha (bottom, right) filters show a larger field of view than the images in the left column. The 
artifacts noticeable around the field stars in the continuum subtracted Br7 and H2 images are residuals of the subtraction caused by a 
slight mismatch of the PSF between the line and continuum images. The arrows towards the northwest of the nebula in the right panels 
mark diffuse emission in the H2, [N ll], and Ho; images associated with the northwest H2 lobe and its northwest extension. The arrow 
towards the south of the nebula in the H2 image marks diffuse emission unrelated to NGC 6881, probably a background galaxy. 
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Fig. 2. WHT LIRIS H2 (red), and NOT ALFOSC Ra (blue) 
and [N ll] (green) composite picture of NGC 6881 overlaid with the 
different morphological features identified in the nebula. 



([Guerrero fc Manchadol [19981 ). The details of the molec- 
ular and ionized material distribution in the central re- 
gions of NGC 6881, how e ver, c ould not be studied prop- 
erly by iGuerrero et alJ (|2QQQ| ). because of the limited 
spatial resolution of their near-IR images and poor con- 
tinuum subtraction at this region. Our new continuum- 
subtracted H2 image reveals two bright knots of H2 emis- 
sion located at the central region of NGC 6881 along its 
minor axis, suggesting a ring of molecular material. The 
relative distribution of ionized and molecular material in 
the central region of NGC 6881, shown in Figure [3l re- 
veals that the H2 emission encompasses the [N ll] emis- 
sion. 

Molecular hydrogen in this region survives in a thin 
layer surrounding the ionized torus-like structure, most 
likely shielded by dense material from the ionizing flux 
of the central star. 

3.1.2. The Highly Collimated Bipolar Lobes 

The two pairs of highly collimated bipolar lobes are 
mainly detected in the Br7, [N ll], and Ha images, while 
the H2 emission from these lobes is weak. The northern 
wall of the northwest lobe is brighter than its southern 
wall, and the opposite applies to the southeast lobe, thus 
showing the point-symmetric brightness distribution typ- 
ical of other bipolar PNe (e.g., K4-55, ■ Guerrero et all 
[1996 ). 

3.1.3. The Hourglass H2 Bipolar Lobes 

Contrary to the highly collimated bipolar lobes, the 
open hourglass bipolar lobes of NGC 6881 are dominated 
by H2 emission. 

They show a narrow waist coincident with NGC 6881 
central regions and clear limb-brightening, indicating 
that they are formed by a thin layer of material. 

The southeast lobe, extending ~12'', is smaller than 
the northwest lobe, with a size up to ~20''. There are 
notable differences between the tips of the southeast and 




Fig. 3. — Expanded view of the central region of NGC 6881 
in the HST WFPC2 [N ll] image. The image is overlaid with H2 
contours extracted from the WHT LIRIS continuum-subtracted H2 
image. 



northwest lobes. The tip of the southeast lobe is sharp 
and can be associated with the southeast loop-like struc- 
ture detected in Ha and very prominently in [N ll] (Fig- 
ure [4j). The westernmost emission of the [N ll] loop is 
coincident with the southwestern elbow of the H2 lobe, 
and the [N ll] loop structure traces closely the edge of the 
H2 lobe. The northwest lobe does not have a sharp edge 
nor does it displays a loop-like structure, but only a knot 
and a thin filament visible in the [N ll] image located at 
its Westernmost tip (Figs. [1] and [3]). 

3.1.4. The Northwest Bipolar Lobe Extension 

The H2 emission of the northwest hourglass lobe dis- 
plays an extension that casts its emission outwards up to 
~50'' from the center of NGC 6881 along its symmetry 
axis. This emission is mainly detected in H2, but there 
are also hints of [N ll] and Ha emission (Fig. [2]). The 
H2 surface brightness of this feature is not homogeneous, 
but is distributed on a series of three bright bands in- 
terspersed between bands of diminished brightness. In- 
terestingly, the bright bands can also be recognized in 
the Ha image. The northwesternmost tip of these bands 
shows a remarkable wedge-shaped morphology. 

3.1.5. Large- Scale Emission 

The Ha and [N ll] images of NGC 6881 show large- 
scale diffuse emission, with the nebula laying on a 
broad arc of patchy emission that crosses the FOV 
from the northeast to the south (Fig. [1]). An exami- 
nation of an Ha image from the INT/WFC Photomet- 
ric Ha Survey (IP HAS) of the northern Galactic plane 
(|Drew et al.l l2QQ5l ) of this region reveals a complex sys- 
tem of filaments. NGC 6881 (/=74?5520, 6=-h02?1137) 
is projected along the local spiral arm onto the intri- 
cate network of Ha filaments near the proximity of 7 
Cyg (e.g., Parker. Gull, fc Kirshnerlll979 ). In particular, 
NGC 6881 is projected towards the east of a large cavity 
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TABLE 4 

Line Fluxes JH and K^ Grisms at PA 113° 



Fig. 4. — Expanded view of the southeast loop-hke feature from 
the H2, Ha, and [N 11] composite picture of NGC 6881 (Fig. [2]). 



TABLE 3 

Line Fluxes at P. A. 137° with K^ Grism 



A 

JAL 



Line ID 



Regions 

Central Ionized H^ NW H2 lobe 

Regions lobes lobes extension 

[xlO^^ ergs cm-2 ^-^ A^^ ] 



19446 


Br8 H I 


63 


19549 


He I 


16 


19574 


H2 (1,0) S(3) 


12 


19703 


H2 (8,6) 0(2) 


2.1 


20338 


H2 (1,0) S(2) 


4.2 


20377 


? 


4.7 


20587 


He I 


47 


20732 


H2 (2,1) S(3) 


1.2 


20763 


7 


0.9 


21126 


He I 


5.1 


21138 


He I 


2.2 


21218 


H2 (1,0) S(l) 


12 


21542 


H2 (2,1) S(2) 


0.4 


21614 


He I ? 


2.9 


21658 


Br7 H I (Brg) 


110 


21793 


Br I 


1.1 


21887 


He II 


17 


21985 


? 


2.1 


22235 


H2 (1,0) S(0) 


3.3 


22477 


H2 (2,1) S(l) 


1.6 


22872 


Br I ? 


7.1 




A 



Regions 



Line ID 



Ionized Lobes H2 Lobes 

^15 



[xlO-*^^ ergs cm 



10936 
11621 
11666 
11892 
11970 
12566 
12784 
12817 
13202 
15335 
15549 
15693 
15875 
15995 
16102 
16401 
16432 
16763 
16801 
16997 
17356 
17478 
19570 
19703 

■ 20338 

■ 20580 
20656 
20732 

■21133 
21218 
21658 
21892 
21986 
22014 
22235 
22477 
22870 



Pa6 H I 

C I ? + Fe I ? 

H2 (3,1) S(4), He I ? 

H2 (2,0) S(0) 

He I 

[Fe II] ?, + Fe II ? 

He I 

Pa 5 (PaB) H i 



H I 4-15 



H I 4-13 
H I 4-12 



Br 18 H I 
Br 16 H I 
Br 15 H I, 
Br 14 H I 

[Fe II] 
Br 13 H I, 
Br 12 H I, 
[Fe II] 
? 

Br 11 H I, H I 4-11, Fe i 

He I 

Br 10 H I 

H2 (1,0) S(7) 

H2 (1,0) S (3) - 

H2 (8,6) O (2) 

H2 (1,0) S(2) 

He I 

H2 (3,2) S(5) 

H2 (2,1) S(3) 

He I 

H2 (1,0) S (1) 

Br 7 (Brg) H I 

He II 

? 

H2 (3,2) S(3) 
H2 (1,0) S(0) 
H2 (2,1) S(l) 
H2 (3,2) S(2) 



[Fe 



160 
21 
4.5 
21 
4.0 
33 
13 
190 
13 
2.3 
1.2 
1.2 
1.1 
1.0 
2.5 
2.6 
15 
1.7 
3.9 
2.5 
5.5 
1.9 
17 
1.1 
3.0 
4.6 
0.5 
0.8 
0.7 
6.8 
7.4 
0.9 
0.7 
1.0 
2.3 
1.2 
1.0 



2 
11 



10 
61 

7.7 



1.6 
11 
0.9 

1.0 

4.7 
1.7 
1.8 



of size 20^x12' oriented along the north-south direction 
that forms part of the H ll region Sh 2-109 (/ = 79?48, h 
= +00?15). 

3.2. Spectroscopic Analysis 

The long-slit spectra of NGC 6881 detect continuum 
emission from the innermost regions and extended line 
emission all through the nebula. A preliminary inspec- 
tion of the two-dimensional spectra reveals the variation 
at different locations of the nebula of important line ra- 
tios (e.g. H2 1—0 S(l)/Br7, He l/He ll), in agreement 
with the different spatial distribution of emission shown 



by the Br7 and H2 images. Guided by these images, we 
have extracted spectra from four individual regions: the 
central region, the ionized lobes, the H2-dominated hour- 
glass lobes, and the northwest H2 lobe extension. The K 
spectra of these regions are shown in Figure \5\ and the 
JHK spectra of the ionized lobes in Figure[6l Within our 
limited spectral resolution, the radial velocity derived for 
the different regions is similar. 

The measured line intensities are listed in Table 3 and 
-4. The spectra and line ratios confirm that the central 
region is Br7-dominated, the ionized lobes show H2/Br7 
~ 1, and the outermost regions are H2-dominated. He ll 
emission is confined to the central region, while the out- 
ermost regions show a wealth of H2 lines. We note that 
the spectrum from the ionized lobes is most likely con- 
taminated by emission from the H2-dominated hourglass 
lobes. 



3.3. Hn Excitation 



The molecule of H2 can be excited by shocks or by 
UV fluorescence. The 1-0 S(l)/2-l S(l) line ratio is 
traditionally used to diagnose the H2 excitation mech- 
anism, with a line ratio ~2 implying UV fluorescence 
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Fig. 5. — Kb spectra along P. A. = 137° of the central region, ionized lobes, H2-dominated bipolar lobes, and northwest extension of 
the northwest H2 bipolar lobe of NGC6881. Line identifications are overlaid on the spectra. 
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lobes of NGC 6881. Line identifications are overlaid on the spectra. 



(|Black fc van Dishoeck|[T987), and a line rat io ~10 im- 
plying shock excitation ([Burton et alJ fToSOl ). The 1- 
S(l)/ 2-1 S(l) line ratio in the different regions of 
NGC 6881 ranges from 5.5 up to 7.5, thus indicating that 
shock excitation can be the most likely dominant excita- 
tion mechanism. 

In a case of pure collisional excitation, the population 
of the energy levels of the H2 molecule are described by 
a single-temperature Boltzmann distribution. To verify 
this statement, the excitation diagrams for several re- 
gions of NGC 6881 (Figure [7j) have been calculated us- 
ing a Boltzmann distribution of states where the column 
densities are related to the ^^ = 1, J = 3 state by 



gsNiu, J) 



exp 



E{iy,J)-E{l,3) 



(1) 



gj and gs being the statistical weights. This equation 



can be rearranged in terms of the flux 

F{u\ JO^l, 05(1)^1, 05(1)^3 



F(l, 3)z/Azy,AJ^z.',J'^r.'',J''^J 



exp 



E{iy,J)-E{l,3) 

K/J- ex 



(2)' 

where Ayf^j>^y>> j,, represent the transition probability, 
and Fiv' ^ J') is the observed flux at frequency iyAu,Aj- 

The temperature can be inferred from the slope of the 
line fltted to the data points in Fig. [71 In NGC 6881, the 
value of the vibrational excitation temperature, T^xiy)^ 
is similar to this of the rotational excitation temperature, 
Tex{J), and thus it can be concluded that shocks are the 
dominant excitation mechanism of H2. The excitation 
temperature is 2100 K for the central region, 2200±100 K 
for the ionized and H2 lobes, and 2700±100 K for the 
extension of the NW H2 lobe. There is thus, a trend for 
the temperature to increase as we move farther from the 
central regions of NGC 6881. 

4. DISCUSSION 

4.1. Ionized and Molecular Gas in NGC 6881 

The II2 emission of bipolar PNe is typically found in 
their equatorial rings and on the walls of their bipolar 
lobes, outlining closely the distribution of ionized mate- 
rial shown in Ha o r [N ll] images (Kas tner et al.lfl996l : 
iGuerreroet all [2000 ). In NGC 6881, however, the differ- 
ent spatial distributions of ionized material and molecu- 
lar hydrogen trace distinct nebular structures. A search 
in the literature shows that there are very few cases of 
bipolar PNe and proto-PNe in which the H2 and Ha 
or [N 11] morphologies are different: CRL2688, M2-9, 
NGC 2440, NGC 7027, J 900, and Hb 12. From this short 
list, NGC 2440 has to be discarded as recent HST images 
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(P roposal ID. 1 1 090) have shown that the H2 arc seen 
by iLatter et alJ (|l995l ) follows the outermost Ha emis- 
sion. In CRL2688 (the Egg Nebula), the H2 emission 
is distributed along the equatorial plane an d at the tips 
of th e bipolar lob es seen in scattered light (Latter et alj 
119931 : ISahai et all 1 1998 ). The H2 bipolar lobes of M 2-9 
envelop the bipolar lobes seen in emission lines of ion- 
ized species such as Fe ll (jSmith et al.ll2005i) . The four- 
lobed shell of H2 in NGC 7072 surrounds its elliptical 
ionized core, revealing the location of a photodissocia- 
tion region ( PDR) that is exci t ed by the absorpti on of 
UV photons (jHora et all Il999l : ILatter et al]l2000f ). In 
J 900, the H2 emission shows a linear struc ture extend- 
ing ^40^^ that ends in an arc-like structure (|Shupe et al.l 
1 19951 : ILatter et al.l Il995l ) . This feature has been sug- 
gested to be caused by UV excitation of stellar radia- 
tion that leaks out thr ough holes in the nebular enve- 
lope (jShupe et al.lll995l ). Finally, the inner ionized bipo- 
lar lobes of Hb 12 are surrounded by an H2 eye-shaped 
structure that is interpreted as a bro ad equat orial ring 
and a larger pair of bipolar lobes (Hora fc Latter 199S 
IWelch et al.l[l999l : iKwok fc Hsiall200l ^ ~^ 

Therefore, among the "exclusive" group of PNe with 
different distributions of H2 and ionized gas, M2-9, 
Hbl2, and NGC 6881 have ionized bipolar lobes that 
are enclosed within II2 bipolar lobes. We note, how- 
ever, that Hb 12 i s a notable case of UV- excited fluores- 
cent H2 emission (|Luhman fc Riekdl 19961 : iHora fc Latter 
[1996) and the same is suspected for M 2-9 (jSmith et al. 
[2005), while NGC 6881 is representative of collisionally 
excited H2 emission. The presence of multiple pairs of 
bipolar lobes in PNe is further discussed in the following 
section. 

4.2. Multiple Bipolar Ejections in NGC 6881 

The occurrence of multiple ejections of material is com- 
mon among PNe, particularly among type I multiple- 
shell PNe that show faint, giant haloes (|Chu et al.l 
Il987l ). The haloes, in terpreted as the relic of late 
thermal-pulse episodes (Stanghellini fc Pasquali' 'l995j; 
[Hajian et al. 1997), are mainly associated with elliptical 
or round PNe, but are rare among bipolar PNe. Recur- 
rent bipolar ejections resulting in multiple independent 
pairs of bipolar lobes have been reported for only a hand- 
ful of PNe. The different pairs of bipolar lobes may have 
simil ar morphology and o rientation, e.g., M2-9 (Balick 
[1999I : ISchwarz et al.iri997l ). but in most cases they show 
notable differences, either affecting their morphology and 
degree of collimation, e.g., M z3 ([Guerrero et al.l [20041 : 
ISantande r-Garcia et al.l l2004l ). the orientation of their 
symmetry axes, e.g., M2-46 |Manchado et al. 1996), or 
its ionization degree, e.g., Hb 12 ( Welch et al .11 19991 ). In 
NGC 6881, we find two pairs of highly collimated ion- 
ized lobes oriented along slightly different directions and 
a pair of H2 bipolar lobes with an open hourglass mor- 
phology. 

The two p airs of ionized l o bes o f NGC 6881, first 
reported by iManchado et al.l (I1996D. hav e been de- 
scribed in detail b y iGuerrero fc Manchadol ((1998) and 
iKwok fc Sul (|2005l ). Precession is evidenced by the dif- 
ferent symmetry axes of each pair of ionized lobes and 
from the alignment between the equatorial ring and the 
youngest pair of bipolar lobes. The sharp morphology, 
similar kinematical age, and different orientation of the 
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Fig. 7. — Excitation diagrams for the central region (top), 
ionized lobes (center), and H2-dominated lobes and northwest ex- 
tension (bottom) of NGC 6881. The diagrams show the upper 
state vibration-rotation populations relative to the 1^ = 1, J = 3 
level plotted against the upper state energy. The linear fits to 
the populations of different vibrational levels find rotational exci- 
tation temperatures of 2000 K, 2100 K, and 2700 K for the central 
region, ionized lobes, and H2-dominated lobes and northwest ex- 
tension, respectively. 



two pairs of bipolar lobes point to a sudden collimated 
ejection of material that has carved out bipolar cavi- 
ties into the surrounding medium. The morphology of 
the ionized bipolar lobes of NGC 6881 is very reminis- 
cent of the proto-PN CRL618 in whi ch hi^h- velocity 
jets moving along different directions (^Cox et alJ l2QQ3l ) 
have produced bipolar caviti es in the nebular envelope 
(|Trammell fc Goodrich! [20021 ). The similarities between 
the two bipolar nebulae extend to the coaxial rings of 
the bipolar cavities and to the bright, compact knots at 
their tips. The collimated outflows of CRL618 might 
well evolved in the future into bipolar lobes similar to 
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those of NGC 6881. 

The hmb-brightened morphology of the H2 lobes indi- 
cates that the material is mainly confined in their walls, 
thus suggesting that the bipolar lobes carry enough mo- 
mentum to sweep up and press the circumstellar medium 
into a thin sheet. Indeed, the prevalent shock excita- 
tion of H2 requires a source of momentum. It is unlikely 
that a current fast (~1000 km s~^) wind from the cen- 
tral star would provide this momentum, as such a stel- 
lar wind would be trapped by the inner ionized lobes. 
Our spectra do not provide information on the expan- 
sion velocity of the H2 lobes, but an upper limit can 
be derived assuming that the [N ll] emission from the 
loop-like structure at the edge of the southeast lobe is 
produced by a shock. For a planar shock to be able to 
produce significant [N ll] A6584 A emission, but low or 
minim al [O ill] A5007 A emission (|Guerrero fc Manchadol 
I1998D, the shock velocity must no t exce ed ~60 km s~^ 
(|Hartigan, Ravmond fc Hartmamil Il987l ). Such high 
speed makes unlikely the origin of the momentum in 
the thermal pressure of hot gas within the bipolar lobes. 
However, the momentum of the II2 bipolar lobes must 
have been provided at the ejection time, either in a sud- 
den bipolar ejection or through a stellar wind intrinsi- 
cally bipolar or isotropic but collimated by an equatorial 
structure. The oblique shock of a bipolar lobe expand- 
ing at ~60 km s~^ would be able to shock excite the H2 
molecules along its edge, while giving rise to [N ll] and 
Ha emission at its tip. Such ballistic ejection may also 
explain the increasing H2 temperature excitation with 
the distance to the nebula center. 

The very different collimation degree of the ionized and 
H2 lobes of NGC 6881 suggests that they formed dur- 
ing two different ejection events with substantial changes 
in the collimation conditions or even in the collimation 
mechanisms between the two episodes. The II2 lobes 
formed very likely by the action of an episodic wind pro- 
ducing an hourglass-shaped ejection, while high velocity 
jets ejected along close directions produced the two pairs 
of ionized lobes. An additional period of isotropic heavy 
mass loss resulting in a dense envelope is revealed by 
concentric arcs detect ed around the innermost regions of 
NGC 6881 (,Corradietl[l.ll2004D . 

4.3. NGC 6881 and the Surrounding Interstellar 
Medium 

The morphology of the southeast and northwest II2- 
dominated bipolar lobes are discrepant, with the former 
being smaller than the latter and showing a sharp edge 
at its southeast tip. This sharp edge is spatially coin- 
cident with a [N 11] loop-like feature (Fig. [4j), a struc- 
ture th at has been interpreted as a lo w velocity precess- 
ing jet ([Guerrero fc Manchadol fl998l ). but that, in view 
of the perfect spatial coincidence with the southeast tip 
of the H2 lobe, seems more likely to mark its termina- 
tion. The abrupt border of the southeast H2 lobe and 
the [N 11] loop-like structure may unveil the interaction 
of this bipolar lobe lobe with dense material. 

As we noted in §3.1, NGC 6881 is embedded within 
diffuse emission that forms part of Sh 2-109, an inho- 
mogeneous H ll region located at a distan ce of 1.4±Q.4 
kpc a nd with a diameter of 400±100 pc (|Fich fc Blitd 
Il984f ). In view of the hint of interaction of the south- 
east H2 lobe of NGC 6881 with dense material, result- 



ing in the notably different morphologies of the south- 
east and northwest lobes, it is worthwhile to investi- 
gate the physical association of NGC 6881 with Sh 2-109. 
The di stance to NGC 6881 has been estim ated to be 2.5 
kpc by Cahn. Kaler fc Stanghellinil (11992 ") who used the 
Shklovsky method according to the Daub scheme. In 
their work, Cahn et al. used an angular radius of 2% 
which seems inadequate as it only traces the bright, cen- 
tral regions of NGC 6881. This region has an angular 
radius ~1'.'2, but if we account for the ionized bipolar 
lobes, then an averaged angular radius of &!5 has to be 
used. If we adopt this radius, then the revised distance 
to NGC 6881 using the Cahn et al. method would be 
~1.5 kpc, supporting its physical association with the 
H II region Sh 2-109. 

An inspection of IPHAS Ha images of this region re- 
veals that the diffuse emission in which NGC 6881 is em- 
bedded constitutes the eastern wall of a large cavity-like 
structure. H I 21 cm line observations of the neutral 
gas distribution and kinematics of the region confirm 
that this struct ure is a cavity with an expansion veloc- 
ity ~10 km s~^ (|Cappa et al.lll996l ). The radial velocity 
of this H I cavity, t'LSR ^ +16 km s~^, is very similar 
to the radial velocity of NGC 6881, ^lsr ^ +2 km s~^ 
([Guerr ero fc Manchadol [19981 ) . The projection of the 
nebula onto the H I cavity and their similar radial ve- 
locities suggest that NGC 6881 may be located within 
this cavity, with its H2 -dominated southeast lobe hitting 
onto the dense material at the wall of the cavity, slowing 
it down and producing a sharp edge, while the northwest 
II2 lobe expands into the H I cavity with little opposition. 

5. CONCLUSIONS AND SUMMARY 

New images and near-IR spectra of NGC 6881 confirm 
the different distribution of ionized material and molec- 
ular hydrogen within this PN. The H2 images have re- 
solved its central region which shows the spatial distribu- 
tion expected for a ring of molecular material surround- 
ing an innermost ionized ring. The H2 bipolar lobes share 
the same orientation than the ionized lobes, but they are 
less collimated and display a distinct hourglass morphol- 
ogy. The H2 emission is predominantly shock excited, 
both in the central ring and in the H2 bipolar lobes. 

We infer a complex history formation for NGC 6881, in- 
volving important changes in the collimation mechanism 
of the bipolar lobes in the latest evolutionary stages. The 
H2 bipolar lobes are consistent with the action of a bipo- 
lar wind. The final evolution of the PN implied the ejec- 
tion of episodic fast collimated outfiows with changing 
directions that interacted with the dense circumstellar 
envelope ejected in a previous phase. 

The southeast II2 lobe of NGC 6881 is less extended 
than the northwest lobe and shows a sharp edge, imply- 
ing the interaction of NGC 6881 with an inhomogeneous 
interstellar medium. We find very likely that NGC 6881 
is within a cavity in the midst of the H ll region Sh 2-109. 
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